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Abstract. The effect of confinement on the phase behaviour and structure of fluids made of biaxial hard 
particles (cuboids) is examined theoretically by means of Onsager second-order virial theory in the limit 
where the long particle axes are frozen in a mutually parallel configuration. Confinement is induced by 
two parallel planar hard walls (slit-pore geometry), with particle long axes perpendicular to the walls 
(perfect homeotropic anchoring). In bulk, a continuous nematic-to-smectic transition takes place, while 
shape anisotropy in the (rectangular) particle cross section induces biaxial ordering. As a consequence, 
four bulk phases, uniaxial and biaxial nematic and smectic phases, can be stabilised as the cross-sectional 
aspect ratio is varied. On confining the fluid, the nematic-to-smectic transition is suppressed, and either 
uniaxial or biaxial phases, separated by a continuous trasition, can be present. Smectic ordering develops 
continuously from the walls for increasing particle concentration (in agreement with the supression of 
nematic-smectic second order transition at confinement), but first-order layering transitions, involving 
structures with n and n + 1 layers, arise in the confined fiuid at high concentration. Competition between 
layering and uniaxial-biaxial ordering leads to three different types of layering transitions, at which the 
two coexisting structures can be both uniaxial, one uniaxial and another biaxial, or both biaxial. Also, the 
interplay between molecular biaxiality and wall interactions is very subtle: while the hard wall disfavours 
the formation of the biaxial phase, biaxiality is against the layering transitions, as we have shown by 
comparing the confined phase behaviour of cylinders and cuboids. The predictive power of Onsager theory 
is checked and confirmed by performing some calculations based on fundamental- measure theory. 



PACS. XX.XX.XX No PACS code given 

1 Introduction 

The behaviour of confined nematics is a problem of great 
fundamental and practical interest. Nematic materials ex- 
hibit strong responses to even subtle external or surface 
fields, which are used in various technological applications. 
. In general, confinement into a narrow planar slit pore 
' induces competing capillary and layering transitions at 
low temperature or high concentration. Recent theoretical 
works have shown that these phenomena lead to complex 
phase behaviour in confined lamellar or smectic phases in 
three [T115J|3] and two @] dimensions. Classical techniques 
[5l[6l[7] and novel atomic-force microscopy have been used 
to study capillary nematization [5112] and the presmectic 
regime |10| . but the strong smectization regime, where lay- 
ering and commensuration effects are crucial, is yet to be 
experimentally studied. 
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Another fascinating topic in the field of liquid crystals 
is the search for biaxial nematic phases ^llj (for a recent 
review, see Berardi et al. [H]). The recent discovery of bi- 
axial nematics in bent-core molecules P^IT^ has renewed 
the interest in biaxiality and in the question of what the 
minimum molecular interaction requirements are for bulk 
biaxial nematic stability. Several factors, such as hydrogen 
bonding, influence the stability of the biaxial ordering |15) . 
The onset of biaxiality and the development of long-range 
bulk biaxial order may also be greatly affected by the pres- 
ence of a surface, a situation where nontrivial coupling be- 
tween biaxiality and surface interactions may occur. Sur- 
faces may induce biaxial phases in one-component fluids 
[161I171IT5] and in mixtures [19 . In addition, confinement 
may induce smectic-like stratification of the fiuid, which 
may be coupled to biaxial ordering. Therefore, the inter- 
play between biaxiality and the capillary nematic-smectic 
transition, and/or layering transitions deep in the smectic 
regime, is expected to yield interesting physics, and is the 
topic of the present paper. 
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Fig. 1. Schematic of particle shape and lengths, and geometry 
used in the confined fiuid. 

Many studies have been devoted to investigate how the 
bulk phenomenology of liquid crystals is affected by strong 
confinement and surface interactions f20l l21l[22] , which may 
promote or discourage layering and biaxiality inside the 
pore. Up to now microscopic models have been scarcely 
used to address this issue. Elastic theory has been em- 
ployed to discuss the anchoring energy in a biaxial nematic 
phase (see e.g. [53]). The phenomenological Landau-de 
Gennes theory has been extensively used to study uni- 
axial and biaxial nematics in films under various condi- 
tions pil[^[^[?7l[^[^[3ni[ni[3^|Ml|Ml[3^[HH] . and also 
specifically in connection with the biaxial transition tak- 
ing place in a thin nematic layer subject to antagonistic 
easy axes in the two surfaces [37l[38l[39ll40] (the so-called 
hybrid cell) . Experimental verification of this phenomenon 
now exists f5T|. The presmectic film [32] and fully devel- 
oped confined smectics [43] were also considered. Mod- 
els dealing only with attractive particle interactions (of 
the Maier-Saupe type) have also been used for confined 
nematic liquid crystals 44,45,46,47j. Theories on lattice 
models [HJHS , 49 and computer simulations on lattice [501 
[5T|[5^[53] and continuous [Ml[55|[H()l[H71[HRl[5!?l[mTl|M] mod- 
els are also abundant. In all these studies the issue of bi- 
axility is implicitely or explicitely considered. 

Density-functional theory (DFT) allows us to study 
the properties of inhomogeneous fluids from molecular in- 
teractions. DFT has been used to study confined nemat- 
ics from a more microscopic perspective 62,63,64,65,66, 
[ST]. A number of studies based on DFT have focused on 
confined nematic phases between symmetric [68ll69ll70l[64] 
and asymmetric [5T] substrates, and also confined smectic 
phases [51I7T] have been studied. All of these DFT studies 
considered only uniaxial, but not biaxial, particles. Since 
DFT gives a direct link between particle shape and phase 
biaxility, it is a very powerful technique to investigate con- 
fined biaxial nematics and will be the tool used in the 
present work. 

All mesogenic molecules are biaxial to some extent. 
Therefore, microscopic studies of idealised biaxial particle 
models in strong confinement and subject to particular 



surface interactions seem pertinent. Following Vanakaras 
et al. [72] , we have recently investigated the bulk phase be- 
haviour of a simplified athermal model for biaxial particle, 
consisting of hard cuboids of lengths a2 < (Ji < L, with 
frozen orientation in the long molecular axes {z axis, see 
Fig.jl]). We used a sophisticated free-energy density func- 
tional based on the fundamental-measure concepts |73) . 
The shape of these particles effectively describes bent-core 
molecules, with the angle of aperture being related to the 
cross-sectional aspect ratio k = cri/cr2 > 1 of the cuboids. 
In this model, there is always a primary nematic director 
z, parallel to the particle long axes; but a second director 
may arise in the xy plane, associated with a nonuniform 
distribution in the azimuthal angle ip. In this case a biax- 
ial phase appears. In line with the findings of Vanakaras 
et al. [7^, based on Onsager theory, we found, in both the 
nematic and smectic regions of phase stability, uniaxial 
and biaxial phases, with a complex bulk phase diagram 
exhibiting continuous transition curves meeting at a four- 
phase point. 

Here we use the same Onsager theory as in Ref. |72) . 
since it is much simpler than our previous theory |73) , but 
predicts essentially the same qualitative behaviour. The 
bulk phase diagram in the plane chemical potential fj, (or 
scaled particle concentration cq) versus aspect ratio of the 
particle cross section k, is shown in Fig. [5] (see later for 
a more comprehensive discussion). The four bulk phases 
identified are: uniaxial nematic (N), biaxial nematic (Nb), 
uniaxial smectic (S), and biaxial smectic (Sb). Depending 
on the value of k, biaxiality may arise in the nematic phase 
(for large aspect ratios) directly from the uniaxial nematic, 
or require as a prerequisite the formation of smectic layers, 
giving rise to a biaxial smectic. 

In this work hard walls will be used to confine the 
fluid, promoting perfect homeotropic anchoring (long par- 
ticle axes perpendicular to the walls. Fig. |T]). The walls 
do not couple to the azimuthal angle ip, so that the sec- 
ond director, whenever there is one, can be taken to be 
aligned in a fixed direction, say the x axis (i.e. both walls 
have identical easy axes; the case with different easy axes 
will be the subject of a future publication). As will be 
seen later, a hard wall turns out to have a disordering ef- 
fect with respect to biaxiality. As a consequence, biaxial 
phases grow in the central region of the slab, not near the 
surfaces, and occur at higher pressure/chemical potential 
than in bulk. However, as common in adsorption problems 
involving hard interactions, the opposite effect occurs for 
the smectic order parameter, which is reinforced near the 
walls, so that strong layering may occur below bulk sat- 
uration. Layering transitions, associated with commensu- 
ration effects between pore width and smectic period, are 
disfavoured by biaxiality. The capillary biaxiality transi- 
tion interacts with the layering transitions in interesting 
ways, giving rise to three types of layering transitions, in- 
volving two biaxial, one uniaxial and another biaxial, or 
both uniaxial, coexisting structures, differing by a single 
layer. Even though the biaxial bulk transition survives to 
confinement, the continuous bulk nematic-smectic tran- 
sition is suppressed under confinement, so that the con- 



S. Varga et al.: Competition between capillarity, layering and biaxiality in a confined liquid crystal 



3 



4.5 



(a) 




1^ 



* 



1 

(b) 


1 


1 










/ 






/ 






/ 






( 


Nb 












1 1 1 


N 

1 1 


1 1 1 



13 14 15 16 17 18 

K 

Fig. 2. Bulk phase diagram of hard cuboids as obtained from 
Onsager theory, (a) Density co (scaled with second virial coef- 
ficient B2) vs. cross-sectional aspect ratio k, and (b) reduced 
chemical potential /i* vs. k. Curves show the boundary be- 
tween the different phases. Dashed curves indicate transitions 
that disappear on confinement. All phase transitions are con- 
tinuous. Labels indicate stable phases. The corresponding bulk 
N-S transitions of the two confined fluids studied are indicated 
by circles. 



transitions are continuous and meet at a four-phase point. 
The derivatives of the transition densities with respect to 
the aspect ratio at different sides of the four-phase point 
are different for the four transition lines. The continu- 
ous nematic-to-smectic transitions, N-S and Nb-Sb, are 
suppresed under confinement (this is indicated by dashed 
curves). The uniaxial (U)-to-biaxial (B) phase transition 
(which, depending on the value of k, is associated with 
either the N-Nb or S-Sb bulk transitions) survives (this is 
indicated by the continuous curves) , and becomes a wavy 
line as the wall separation h (see Fig. [1]) is varied; its 
average location is shifted to higher values of \x due to 
the disordering effect of the walls. However, even though 
the bulk nematic-to-smectic bulk transition disappears on 
confinement, it governs the behaviour of the layering crit- 
ical points, since the latter tend to the bulk nematic-to- 
smectic chemical potential as the slit-pore is made wider. 
In one case considered (k = 13), the relevant bulk transi- 
tion is N-S, while in the other (k — 17) layering transitions 
are governed by the Nb-Sb bulk transition (open circles 
in Fig. [2]). In addition, when the U-B transition is above 
the bulk N-S transition (k < 15.304), layering transitions 
intersect with the U-B transition curve, whereas in the 
opposite case the two transitions are independent. 

The outline of the paper is as follows. In Section [5] the 
model system and the theory are presented. The phase di- 
agrams of bulk and confined systems and the equilibrium 
density and order parameter profiles are shown and dis- 
cussed in Section[31 In Section|3]we make some concluding 
remarks. Finally, we present a bifurcation analysis for the 
N-S, N-Nb and U-B phase transitions in the Appendix, 
and an explanation for the emergence of biaxial ordering 
near the S-Ss transition. 



2 Theory 

We study the effect of confinement on the phase behaviour 
of a fluid of biaxial hard particles. The appropriate ther- 
modynamic free energy for confined fluids in DFT is the 
grand-canonical free-energy functional per unit area i7[p]/^, 
which in our case is given by 



n{p\ F{p\ 



A 



A 



dz / d(pp{z, (fi) [p - Vay.t{z, (p)] , (1) 



fined nematic smoothly transforms into a confined smectic 
phase; here confined smectic phases are identified by their 
propensity to undergo layering transitions at high values 
of chemical potential. Each of these (infinitely large in 
number) first-order layering transitions disappears at a 
critical point. 

Two different scenarios, corresponding to two differ- 
ent choices for aspect ratio, will be presented: k = 13 
and K = 17. To explain these two choices, we again use 
Fig. [2] which represents the bulk phase diagram in the re- 
duced chemical potential ~ l3fi, or scaled particle con- 
centration Co, versus aspect ratio k plane. All bulk phase 



where p(z, ip) is the position- and orientation- dependent 
local number density, and Vc-^tiz, (p) is the external poten- 
tial. In Eqn. (1) the Helmholtz free-energy functional is 
the sum of ideal and excess contributions, F[p] = Fi^[p] -\- 
-Pcx[p], and can be written as 

^ I dzj dMz.v) [\og{p[z,^)A^) - 1] , 



A 



-i j dzi J dip I j dz2 J dp2p{zi,pi) 



'Xp{z2,V2)f{zi2,(pi2), 



(2) 
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where A is the thermal wavelength and f{zi2,(pi2) is the 
integrated Mayer function over the x — y plane, with zi2 — 
z\ — Z2 and ipvi = '^\—'^2- We set A = 1 since this does not 
affect the fluid phase behaviour. Note that, while the ideal 
contribution is exact, the excess part is approximated at 
the level of second-order virial theory. Particles are con- 
fined between two parallel hard walls, the effect of which 
is modelled by an external potential. 



oo, z < 0, 

/?T/ext(^,</?) - { 0, 0<Z<h, 



(3) 



oo. 



> h. 



where h + L is the distance between the walls, and the 
normal to the wall is along the z direction. This wall- 
particle interaction confines the centre of mass of the par- 
ticles to stay in the finite interval < z < h. Since the 
cross section of the cuboid is a rectangle, it is sufficient 
to consider the interval < ip < n. Due to the exclusion 
interaction between the wall and the particles, we have 
p(z, ip) = for z < and z > h. The use of frozen orienta- 
tions for the long particle axes substantially simplifies the 
calculations. This approximation may be valid in a very 
dense mesogenic system with perfect homeotropic anchor- 
ing (homeotropic anchoring could be achieved by means of 
very strong external electric/magnetic fields or by special 
surface treatments leading to specific wall-particle inter- 
actions). After all these assumptions, the grand canonical 
free energy functional per unit volume can be simplified 
to give 



f2[p] _ F[p] p r'^ 



V V 



(4) 



where 



p{z) = / dipp{z,ip). (5) 
Jo 

The ideal and excess free-energy contributions become 
= T / / d(pp{z, ip) {log ip)] - 1} , 



Jo 



T7 / "-^^ / dipip[zi,ipi) / dz2 

y Jo Jq Jzi-a(zi) 



X / dip2p{z2,ip2)Ac^{zi2,(pi2), 

Jq 

where the functions a(z) and b{z) are defined by 



a{z) 




Hz) 




(6) 



(7) 



(8) 



Acx{z,ip) is the excluded area between the cross sections 
of two particles. Functional minimization of Eqn. (|4]) with 
respect to the local density at fixed chemical potential 
provides the equilibrium density profile of the fiuid. The 
resulting integral equation for p(z, 1^9) is 

I3p- dz2 / dp2P{z2,'f2)Acx{zi2,pi2) 

— g Jzi-a{zi) Jo j^Q^ 

The inputs to the above equation are the chemical po- 
tential fi and the wall separation h. The excluded area 
between two cuboids is given by 

Acx{z, (fi) — \2a\G2 (1 + I cos(y5|) + {g\ + ctj) I sini^j] 



x0(i-|z|). 



(10) 



We have solved Eqn. © using two methods. One is the 
standard iterative procedure, which was applied in the 
case of wide pores. The other is the Fourier expansion 
method, used for narrow pores. The choice of method was 
dictated by the need to optimise computation time, but 
both give numerically identical results in a wide range of 
pore sizes. 

The properties of the bulk smectic fiuid can also be 
calculated by solving Eqn. © at fixed chemical potential 
and fixed value of h [in the case of the nematic, we simply 
evaluate dl])-® for p{z,ip) = po = N/V, the mean den- 
sity]. For the bulk phases, the a{z) and b{z) functions have 
to be defined as a{z) = b{z) = L, and periodic boundary 
conditions are applied for p{z, (f) along the z direction in 
the case of the smectic phase. The value of h that min- 
imizes the grand potential is identified with the smectic 
period d, while the corresponding density distribution is 
the equilibrium density distribution of the smectic phase. 
In the following, the particle length (long axis) L is used 
as unit length, so that the scaled pore width and the z 
coordinate are h* — h/L and z* = z/L, respectively. In 
the case of density, a more convenient unit of volume is 
the second virial coefficient of the particles in the uniax- 
ial nematic phase, which is given by B2 = LA^^ , where 
A^^ = Tr~^ dip AcxiO,ip). The scaled bulk density will 
be Co = poB2- With this choice the bulk uniaxial nematic- 
smectic transition density is independent of the geometry 
and shape anisotropy of the cross section (see lA.ll) . 

Eqn. ^ furnishes the density profile at given chemical 
potential for a particular value of transverse aspect ratio. 
This density profile can be either uniaxial or biaxial, both 
in bulk and in the confined cases. It may also happen that 
more than one solution exists for a given input. Therefore 
an analysis of the grand potential is necessary to find the 
stable phase and possible phase transitions. In the next 
section we present the bulk and confined phase diagrams, 
as well as density profiles, as obtained from the method 
explained above. 



S. Varga et al.: Competition between capillarity, layering and biaxiality in a confined liquid crystal 



5 



6 
5 
4 

C 3 
2 
1 






1^* = 6.100, / 




1 riQ / 
u.lUo, / 




6.110 / 


(a) 





1 ' r 

1 



3 4 5 6 



5.0 
4.5 
4.0 
3.5 
3.0 
2.5 



0.91 






^*=6.90 








\\ 6.95 // 




\7.00/ 


■(c) 




■ 1 

1 


1 1 1 ■ 1 ■ 1 ■ 1 

2 3 4 5 6 


1 


2Tiz/d 






1^* = 7.00 




/ 6.95 \ 




'/ 6.91 \\ 




6.90^^^^ 


(d) 


■ 1 

1 


1 1 1 ■ 1 ■ 1 ■ 1 

2 3 4 5 6 




2%z/d 



2 3 4 5 
2 7tz/J 

Fig. 3. Density profile and biaxial order parameter in the vicinity of the S-Sb transition ((a) and (b)), and Nb-Sb transition 
((c) and (d)) for hard cuboids. The aspect ratio of the particle cross section is k = 13 in (a) and (b), and k = 17 in (c) and (d). 
The values of the corresponding chemical potentials are shown in the figures. 



3 Results and discussion 
3.1 Bulk behaviour 

The bulk phase diagram of the parallel hard-cuboid fluid 
has already been determined by both Onsager theory ^72j 
and Fundamental-Measure theory (FMT) [75]. The main 
focus of those studies was to determine the effect of the 
cross-sectional aspect ratio on the stability of biaxial or- 
dering against uniaxial ordering. It was found that macro- 
scopic biaxial ordering can be enhanced with increasing 
aspect ratio k. The theories agree qualitatively, in that 
both predict continuous uniaxial nematic-uniaxial smec- 
tic (N-S), uniaxial smectic-biaxial smectic (S-Sb), uniax- 
ial nematic-biaxial nematic (N-Nb), and biaxial nematic- 
biaxial smectic (Nb-Sb) phase transitions, all meeting at 
a four-phase point located at (k4, ^4). However, the tran- 
sition lines and the location of the four-phase point, and 
therefore the onset of Nb stability, are different in the 
two theories; for example, Onsager theory predicts K4 = 
15.304, while K4 = 18.101 with FMT. Nevertheless, it is 
quite surprising that a second-order virial theory such as 
Onsager's works quite well, even in the highly ordered 
smectic phase. Later in this section we show that Onsager 



and FMT theories also produce very similar density pro- 
files in the confined system. 

As a first step, we recalculated the bulk phase dia- 
gram, which is a useful reference to discuss the properties 
of the confined fluid and its connection to bulk behaviour. 
Starting from a high value of /i, and performing minimi- 
sations for decreasing values, the nematic-smectic phase 
boundary for a given aspect ratio k can be obtained when 
the equilibrium density distribution becomes independent 
of z. In addition, we performed N-S and N-Nb bifurca- 
tion analyses (see lA.ll and IA.2] for details), the results of 
which agree perfectly well with those from the minimi- 
sations. Using Co as scaled density, the value of density 
at the boundary, c^^ , is independent of both aspect ratio 
and even shape of the particle cross section. The resulting 
equation for the N-S transition density is 

cT = -[2j,{x*)]-\ (11) 

where tt < x* < 37r/2 is the solution of the trascendental 
equation tana; = x. The phase diagram is depicted in Fig. 
[5] In the calculations that follow, we have chosen two sig- 
nificantly different model systems: k = 13 and 17. As can 
be seen from the figure, these two systems have different 
phase sequences with increasing chemical potential (den- 
sity): N— ^S— >-Sb for K = 13, and N->Nb— 5>Sb for the more 
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Fig. 4. Phase diagram of confined hard cylinders. Phase 
boundaries of the layering transitions in the chemical poten- 
tial /i vs. pore width h plane. The dashed line corresponds to 
the chemical potential of the bulk N-S transition. The layering 
transition in the widest pore examined takes place between the 
structures with 26 and 27 layers. 



anisotropic system {k — 17). In Fig. [3] the local scaled 
density c(z) and biaxial order parameter A(z), defined as 



c(z) 



d(pp{z,(p), 



1 r 

^(^) = -TT / dip cofi2Lp p{z, If), (12) 

pyz) Jo 



are plotted as a function of z along one smectic period. 
The evolution of tliese quantities with fi is quite interest- 
ing. For example, the system with k = 13 undergoes a 
S-Sb transition at fi* = 6.1, but the density profile hardly 
changes with increasing fi, while biaxial order develops 
quickly. A remarkable behaviour of the biaxial ordering is 
that it peaks in the middle of the interstitial region, which 
means that particles are more ordered between neighbor- 
ing layers than inside the layers. A detailed explanation 
for this behaviour, which was already noticed in our pre- 
vious FMT study |73], is provided in Sec. IA.41 Regarding 
the Nb-Sb transition of the k = 17 system [Figs.jHJc) and 
(d)], it takes place at ^* = 6.9, with cq = 3.78 and a biax- 
ial order parameter Aq = 0.886. Both the density and the 
biaxial profiles change substantially on the smectic side. 
As in the previous case, the density and biaxial profiles 
are out-of-phase and, interestingly, the in-layer biaxial or- 
der decreases slightly with increasing chemical potential. 
This is evidently due to the more efficient packing in the 
smectic phase. 
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Fig. 5. Phase diagram of confined hard cylinders, (a) fi*-h* 
phase diagram in the region of the 26-27 layering transition. 
Inset shows the 2-3 layering transition, (b) Reduced chemical 
potential /i* vs. coexistence mean densities cq for the 26-27 
layering transition. Inset shows the 2-3 layering transition. 



3.2 Behaviour under confinement 

3.2.1 Hard cylinders 

We start by discussing a simpler fluid, that of parallel hard 
cylinders. This fluid only undergoes a single transition in 
bulk, the N-S transition, and does not exhibit any biax- 
ial phase (the solid phase is not considered in the present 
study). Due to suppression of long-range smectic fluctu- 
ations, the N-S transition vanishes in the confined case. 
However first-order layering transitions do take place. At 
a layering transition, two smectic-like structures, having 
n and n + 1 layers, coexist at the same pore width, which 
can accommodate slightly compressed or swollen struc- 
tures. The phase diagram is shown in Fig. |4] in the /.i-h 
plane for a wide range of pore widths. Since more and 
more layers can accommodate into the slit pore with in- 
creasing pore width, the number of layering transition is 
infinite in the limit h ^ oo. We can see that the layering 
transition curves are almost straight and vertical in the 
li-h plane. However, the slope decreases slightly with pore 
width. The distance between two consecutive transition 
curves is about 1.35L (with L the length of the cylin- 
ders), and all curves terminate in lower critical points. 
The dependence of the location of these critical points 
on h shows that the layering transitions can be stabilized 
with increasing pore width. In addition, values of the crit- 
ical chemical potentials (or densities) converge to the bulk 
N-S value /iNs with increasing pore width. The shape of 
the layering transition curves is highlighted in Fig.lSJa) by 
showing two extreme cases: the regions around the 2 — 3 
and 26 — 27 layering transitions. The transition curves are 
very steep, but they are convex in very narrow pores and 
concave in wide pores. The dependence of the coexist- 
ing average densities with p is shown in Fig. [5l^b), where 

Co is defined as cq = h^^ /q dzc{z). It is clear that the 
density gap between coexisting layered structures shrinks 
with increasing pore width. This is due to the fact that 
the contribution from wall-particle interactions decreases 
relative to that from particle-particle interactions, i.e. the 
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wall has less effect on the structure in the middle of the 
pore. 




21.0 



Fig. 6. Phase diagram of confined hard cuboids in the ^*-h* 
plane. Layering transitions are indicated by continuous curves, 
while the uniaxial-biaxial phase is indicated by a dashed 
curve. Inset shows the uniaxial-biaxial phase transition for 
narrow pores. U and B denote unia^xial and biaxial phases, 
respectively. Diamonds indicate the intersection between the 
uniaxial-biaxial and layering transitions. The value of the cross- 
sectional aspect ratio is «; = 13. 



3.2.2 Hard cuboids 

In the case of confined hard cuboids again there is no 
N-S transition. However, by contrast with the previous 
case, there is a uniaxial-to-biaxial (U-B) phase transition. 
Therefore only one type of (continuous) ordering transi- 
tion can take place. The phase diagram of confined cuboids 
is shown in Fig. [5] While the layering transition curves 
are obtained by searching for two solutions of the Euler- 
Lagrange equation, Eqn. ©, with equal grand potentials, 
the U-B phase boundary was obtained from the numer- 
ical solution of the set of bifurcation equations derived 
in rOl [see Eqns. (UHl), and The U-B transi- 

tion line ^ub(^) is an oscillatory function of pore width 
that asymptotically tends to the bulk S-Sb transition, 
A'sSis — 6.1, with increasing pore width. This means that 
the wall separation strongly influences the biaxial order- 
ing and, except for some specific intermediate ranges of 
pore width, confinement destabilizes biaxial ordering in 
general. For h > 6.6 the curve /iusCi) interacts with the 
layering-transition structure and becomes discontinuous. 
In this case biaxiality changes the structure of the two co- 
existing phases at a layering transition. With increasing 
fi along a particular layering transition, first the phase 
at left of the transition curve becomes biaxial, and then 



both phases show biaxial order. Therefore three types of 
layering transitions (U-U, U-B or B-B) can emerge in the 
pore. 

The density and order parameter profiles for two coex- 
isting structures of types U-B and B-B at the 5-6 layering 
transition are illustrated in Fig. [T] Due to conimensura- 
tion effects, the density distribution of 5 layers is always 
less peaked than that of 6 layers. In the 5-layer structure, 
the pore width is a bit too spacious for the rods, while for 
the 6-layer structure it is too narrow. This is the case for 
the two coexisting structures at the general n-n + 1 tran- 
sition. Similar to the bulk biaxial smectic phase [see Fig. 
[5], particles are less ordered orientationally at the layers 
than in the interstitial regions. As the wall reduces biaxial 
order and therefore destabilizes the U-B phase transition, 
the biaxial order parameter is lowest at contact with the 
wall. 

The interplay between biaxiality and layering transi- 
tion is revealed in Fig. \E\ where we demonstrate the desta- 
bilising effect of increasing the cross-sectional aspect ratio 
on the biaxial ordering for the 5-6 layering transition. Not 
surprisingly, the uniaxial-biaxial boundary curve moves 
down in chemical potential due to the increased molecular 
biaxiality. But, in addition, there is a reduction in mean 
density gap at the layering transition on increasing k from 
12.7 to 13.3. Also interesting is that the uniaxial phase is 
reentrant in a narrow range of pore widths as /i or cq is 
increased, with a phase sequence U5 —^65 — ^Ug — S-Bg. 
The reentrant phenomenon can be seen most clearly in 
the case k = 12.7. Finally, as can be observed in Fig.|Sl[c), 
the U5-U6 transition disappears by increasing the aspect 
ratio. 

The phase diagram for the k — 17 fluid is shown in Fig. 
ini The U-B phase transition curve /iuB(^) shows damped 
oscillatory behaviour, and in this case converges to the 
bulk N-Nb value with increasing pore width. The desta- 
bilization effect of the walls on the U-B phase transition 
can be seen very clearly in this case. Since the bulk fluid 
now has a large region of Nb stability, and the layering 
transition critical points converge to /iNsSB > /^nnb (see 
inset of Fig. |9]), the U-B boundary does not intersect the 
layering transition curves, which now involve two biaxial 
phases and end in lower critical points. 

Finally, we make some remarks about the reliability of 
Onsager theory for confined systems. It is true that the 
second-order virial approximation seems to be a crude ap- 
proximation for dense non-uniform fluids. To assess the 
validity of Onsager theory, we have made a comparison 
between the predictions of our theory and those of our pre- 
vious FMT approximation [73] formulated for the inhomo- 
geneous fluid. In Fig. [TU]we focus on a particular (15-16) 
layering transition; the density proflles of the two coexist- 
ing structures at the layering transition, as obtained from 
Onsager and FMT theories, and at a given value of ji, are 
shown. In Fig. [TT] the locations of a few layering transi- 
tions in the fi-h phase diagrams are compared, while Fig. 
[12] shows density and biaxial order parameter proflles for a 
given chemical potential. Overall, the agreement between 
the two theories is quite fair. This can be understood if 
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Fig. 7. Density and order-parameter profiles of coexisting structures with 5 and 6 layers. Panels (a) and (b) show profiles at 
the U-B layering transition with /u* = 6.1, while panels (c) and (d) correspond to those at the B-B transition with /it* = 6.3. Us 
and Ue denote uniaxial phases with 5 and 6 layers, respectively, while B5 and Be pertain to biaxial phases. The cross-sectional 
aspect ratio is « = 13. 



we take into account that for high aspect ratios the On- 
sagcr theory usually gives correct phase behaviour. The 
oscillatory structure predicted by Onsager theory is over- 
estimated, while the width of the density peaks is a bit 
wider than it should be. These discrepancies are due to 
the poor treatment of correlations in Onsager theory. The 
location of the layering transitions is shifted in pore width. 
However, Onsager theory can be trusted as far as quali- 
tative behaviours are concerned, and the predicted phase 
diagram is expected to be qualitatively correct. 

4 Conclusions 

We have examined the effect of confinement on the phase 
behaviour of a fluid of parallel biaxial hard cuboids. Using 
Onsager theory, four phases are observed in the bulk limit, 
namely uniaxial nematic, uniaxial smectic, biaxial nematic 
and biaxial smectic. For shape anisotropy with k < 15.304, 
the sequence N— >S— >Sb is observed with increasing chemi- 
cal potential, while the uniaxial nematic phase transforms 
directly into the biaxial nematic phase for k > 15.304, i.e. 
in this case the phase sequence is N^'Nb^-Sb. 

In the fluid confined into a slit pore only two phases 
occur: one has uniaxial structure, while the other is biax- 
ial. The N-S transition cannot survive, but layering tran- 



sitions do exist and uniaxial-biaxial phase transitions do 

occur. Layering transitions are associated to commensu- 
ration effects between a structure with an integer number 
of layers and a finite pore width, while biaxial ordering 
is the result of excluded- volume interactions taking place 
between the cross sections of the biaxial particles. Inter- 
estingly, layering transitions can be destabilized with re- 
spect to an increasing particle biaxiality, characterised by 
a shape parameter k, which means that particles with a 
circular cross sections are the best candidates for observ- 
ing layering transitions. 

Layering transitions between uniaxial-uniaxial, uniaxial- 
biaxial, or biaxial-biaxial, structures can occur at high 
chemical potentials, and reentrant behaviour of the uniax- 
ial phase may occur for increasing chemical potential and 
for some values of pore width. Confinement does not en- 
courage the formation of biaxial order, since the uniaxial- 
biaxial phase boundary moves to higher chemical potential 
(density) for decreasing pore width. Therefore, for real flu- 
ids that can be modelled by this particular hard wall/fluid 
model system, such as colloidal suspensions of anisotropic 
particles, experimental detection of the biaxial nematic 
phase will be more favourable in bulk than in confined 
geometry. 
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Fig. 8. Phase boundaries of layering (continuous curves) and U-B (dashed curves) phase transitions. Upper panels (a) and 
(c): fi*-h* plane. Lower panels (b) and (d): co-h* plane. U5 and Ue denote uniaxial phases with 5 and 6 layers, while B5 and 
Be correspond to biaxial phases. Left panels (a) and (b): k = 12.7; right panels (c) and (d): k = 13.3. 



Our results can be compared to those obtained earlier 
by van Roij et al. P^I17| . who used an Onsager model 
in the Zvifanzig (restricted orientations) approximation on 
hard parallelepipeds of dimensions Lxa x a, with L 3> (t, 
to analyse surface and capillary behaviours. In the Zwanzig 
model particles can only point along the x, y or z direc- 
tions, and the fluid can accordingly be considered as a 
mixture of three species. If we take L = cri, our model can 
be related to the one in Refs. P^I17) , with the important 
difference that, in our case, the species perpendicular to 
the walls is missing, so that the isotropic phase cannot oc- 
cur and there is no capillary nematization/isotropization 
in our model (by contrast, particles can freely rotate in the 
xy plane, and our particle volume is finite, so that layered 
interfacial profiles can be obtained). As a consequence, our 
U-B transition occurs inside the nematic region, whereas 
van Roij et al. obtain it below capillary nematization (i.e. 
in a nematic film adsorbed on the walls), and consequently 
there is no dependence of the U-B transition line with 
pore width (in contrast with our strong oscillatory de- 
pendence, cf. Figs. Eland Our work is a step forward 
in the study of confined fluids of hard particles, in that 
layered smectic phases, not contemplated before, are in- 
cluded, the relationship between biaxiality and layering is 
ascertained, and the dependence of these phenomena with 
particle cross-sectional aspect ratio is analysed. 



Finally, our study shows that Onsager theory can be 
applied for bulk and confined studies even in very ordered 
phases like the uniaxial or the biaxial smectic phases. 
Comparison of the results of Onsager with FMT theories 
indicates that, by including higher-order correlations into 
the theory, only quantitative, rather than qualitative, im- 
provement can be achieved. We think that further studies 
should be performed in order to get a deeper understand- 
ing of the delicate interplay between the surface and biax- 
ial interactions on the stability of layering, capillary and 
orientational transitions. It is not clear how the scenario 
would change if planar, instead of homeotropic, anchoring 
is favoured, or if the orientational entropy is included by 
studying fluids of freely rotating biaxial rods. 

Support from Comunidad Autonoma de Madrid (grants 
S-0505/ESP-0299 and NANOFLUID), Spain, and grants 
FIS2007-65869-C03-01, FIS2008-05865-C02-02 and MO- 
SAICO of the Ministerio de Educacion y Ciencia of Spain 
are acknowledged. 



A Bifurcation analysis 

In this section we provide some analytical results for the 
N-S and N-Nb phase transitions of the bulk system, and 
also for the location of the U-B transition in the con- 
fined fluid. We also give an explanation for the peculiar 
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Fig. 9. Phase diagram in the /i* — h* plane of the U-B transi- 
tion for hard cuboids of aspect ratio n — 17. Upper inset shows 
detail of the same transition for small pore widths. Lower inset 
shows layering transitions ocurring well above the U-B transi- 
tion; horizontal dashed line corresponds to the bulk Nb-Sb- 

behaviour of the biaxial order parameter near the S-Sb 
transition. 



A.l N-S bifurcation 

The bifurcation analysis of the N-S phase transition of 
parallel hard cylinders has been discussed in earlier work 
|74j . Here we only present a short overview of the bifur- 
cation analysis and show the resulting equations for the 
bulk bifurcation density and wavenumber of the present 
biaxial hard particles. The starting point is the free en- 
ergy of the smectic phase, given by Eqns. (H]), for a weak 
periodic perturbation of amplitude e: 

p{z) ^ po{l + e cos qz) , (13) 

where q is the wavenumber. Substitution in Eqns. (jS]) pro- 
vides the free-energy difference AF between smectic and 
nematic phases to lowest order in e: 

PAF 1 



V 



Po + 2plLA[!l^joiqL) 



(14) 



where jo{x) is a spherical Bessel function. At the N-S bi- 
furcation point, 

Pi A p 

Z\F(po,9)=0, ^=0, (15) 



dq 



(16) 



and the transition density results in 

= - [2S2jo(<ZNsir\ 

where tt < q^gL < 3tt/2 is the smectic wave number at 
bifurcation, which can be obtained as a solution of the 
trascendental equation tangNS^ = CfNsL- The dimension- 
less density Cq^ = B2Pq^ is insensitive to the particular 
particle model used. 




Fig. 10. Comparison of the U15-U16 layering transition den- 
sity profiles from Onsager and FMT theories for k = 13 and 
Pp' = Pp + \ogvo = 3.5. (a): Onsager theory, (b): FMT theory. 
The dimensionaless local density is defined as rj{z) — p{z)vo, 
with vo ~ crio"2-L. 



A. 2 N-Nb bifurcation 

Similar to the N-S bifurcation analysis, here we determine 
the free-energy cost associated with an infinitesimal bi- 
axial ordering in the uniform nematic phase. The density 
distribution 



p{lp) ^ — {1 + ecos2(p) , 

TT 



/3AF 
V 



Po 



(17) 



(18) 



with = 27r ^ cos2ipAc^{0, ip). Again the condi- 

tion AF = provides a N-Nb bifurcation density, given 
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Fig. 11. Comparison of the locations of a few layering 
transitions in the n'-h phase diagram for uniaxial confined 
structures. Continuous curves: Onsager theory. Dashed curves: 
FMT. K = 13. A redefined chemical potential n' — /i* + logwo 
was used. 



by 



2LAiV 



(19) 



A. 3 U-B bifurcation in confined fluid 



Unlike the bulk nematic phase, the confined system is 
nonuniform, and the uniaxial local density profile has to 
be determined from Eqn. ^ using a value for the chemi- 
cal potential corresponding to the bifurcation point, which 
is not known a priori. Using Eqn. it is easy to prove, 
by integrating out the ip dependence, that the uniaxial 
density distribution is the solution of the equation 



p{z) = exp 



-a(z) 



(20) 



The biaxial perturbation is again proportional to the func- 
tion cos 2(/3, but its amplitude must depend on the position 
in the pore. Therefore we use the following ansatz for the 
perturbation: 



[I + eF{z) cos2(p] 



(21) 



Substitution into Eqns. ([5]) provides a free-energy differ- 
ence between the uniaxial and the biaxial confined struc- 
tures, which is quadratic in e. The coefficient must be zero 
at the bifurcation point, which provides the equation 



+- 



h rzi'\-b{zi) 

dzif{zi) / dz2f{z2) = 0, (22) 

J zi -a(zi) 




Fig. 12. Comparison of density and biaxial order parameter 
profiles from Onsager and FMT theories for k = 17 and p! — 
4.5 (see caption of Fig. llUp . (a) Density profiles ri{z) = p{z)vo. 
(b) Biaxial order parameter profile A{z). Continuous curves: 
Onsager theory. Dashed curves: FMT. 



where f{z) = F{z)p{z). This equation is still not suitable 
for the determination of f{z), which requires minimizing 
the perturbed free energy with respect to f{z). This cor- 
responds to the functional derivative of with respect 
to /(z) being zero, i.e. 



1 pz+b{z) 

fiz) = -t:A[1}p{z) / dzi/(zi). 



(23) 



The simultaneous solution of Eqns. (HO]), (EH) and (1^ 
gives the density profile p{z), chemical potential p and 
perturbation function /(z) at the uniaxial-biaxial phase 
transition. From the set of U-B bifurcation equations one 
can get analytical results in the case < h < L. This 
interval corresponds to a very narrow pore, where only 
one monolayer of particles can fit in the pore. The upper 
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and lower bounds of the integral in (pn|) become and h, 
respectively, i.e. the integral will give the same value for 
any position in the pore. This means that the local density 
is constant in such a narrow pore. Its value depends on the 
chemical potential and is given by 



p{z) = po = e 



13 p - Al°JhpQ 



(24) 



Similar to the local density, ([25]) predicts a constant value 
f{z) — /o and a density at bifurcation given by 



hA. 



(25) 



Note that ([^5)1 gives the same bifurcation equation as 
/iUB can be obtained from using (^5]) . 



function p{z±L) in general, because we have assumed that 
the smectic phase is highly ordered and thus that the first 
derivative of the density profile p'{z) can reach high values 
for some values of z. Making all these approximations, 
we find the following diferential equation for the order 
parameter: 

^'(^) = IPi^ + - Pi^ ~ L)] . (31) 



The solution to this equation is 



A{z) = Z\(0) exp [If (z) - If (0)], 



tf'(z) 



.(1) 



z+L 



p{z')dz' . 



(32) 



(33) 



A. 4 Smectic phase biaxiality 

This section is devoted to showing how the biaxial order- 
ing of particles develops near the S-Sb transition for the 
particular case where the S phase is highly ordered. We 
explicitly show the peculiar behaviour of the order pa- 
rameter profile as a periodic function peaked around the 
midpoint between adjacent smectic layers. 

Substitution of Eqn. (I^Tl) into the definition of the bi- 
axial order parameter, i.e. 

A{z) — p{z)^^ / p{z, if) cos 2ipdip (26) 
Jo 

provides the result A{z) — F{z)/2. Now insertion of the 
definition f{z) = F{z)p{z) = 2A{z)p{z) into Eqn. ^ 
allows us to find 

A{z) = -i^W 1^^^^ dzip(zi)Z\(zi), (27) 

which can be viewed as an integral equation for the biaxial 
order parameter profile A{z) near the S-Sb transition. The 
first derivative of Eqn. ((27)) with respect to z gives 



A'{z) 



[p{z + L)A(z + L)~ p{z ~ L)A{z ~ L)] 



(28) 



Due to the periodicity of the order parameter profile, the 
function A{z ± L) is periodic with period d: 

A{z±L) = A[z±{L-d)]. (29) 

Its Taylor expansion about z, up to first order, gives 

A{z±L)^ A{z)±A'{z){L-d). (30) 

For the particular case of a highly ordered smectic phase 
(with d ^ L), we can make the approximation A(z±L) ss 
A{z), where use has been made of the fact that A'{z) has 
the same order of magnitude as A{z) near the bifurcation 
point. Note that this approximation is not justified for the 



Assuming that the periodic density profile p{z) is peaked 
at z = Zfe = fed (fc G Z), it is easy to show that the 
maxima of the function tf'(z) given by (j33p are located 
at z = Zk -\- d/2, resulting in an order-parameter profile 
A(z) in antiphase with respect to the density profile p{z) 
[see Eqns. ([32]) and ((33|) ]. The present result is valid for 
any hard non-local interaction with a range of twice the 
particle length. 



References 

1. T. Geisinger, M. MuUer, and K. Binder, J. Chem. Phys. 
Ill, 5241 (1999). 

2. V. Babin, A. Ciach, and M. Tasinkevych, J. Chem. Phys. 
114, 9585 (2001). 

3. D. de las Heras, E. Velasco and L. Mederos, Phys. Rev. 
Lett. 94, 017801 (2005). 

4. Y. Martmez-Raton, Phys. Rev. E 75, 051708 (2007). 

5. H. Yokoyama, J. Chem. Soc. Farad. T. 84, 1023 (1988). 

6. M. M. Wittebrood, D. H. Luijendijk, S. Stallinga, Th. Ras- 
ing and I. Musevic, Phys. Rev. E 54, 5232 (1996). 

7. K. Kocevar, I. Musevic, Phys. Rev. E 64, 051711 (2001). 

8. K. Kocevar, A. Borstnik, I. Musevic and S. Zumer, Phys. 
Rev. Lett. 86, 5914 (2001). 

9. A. Borstnik Bracic, K. Kocevar, I. Musevic and S. Zumer, 
Phys. Rev. E 68, 011708 (2003). 

10. L. Moreau, P. Richetti and P. Barois, Phys. Rev. Lett. 73, 
3556 (1994). 

11. G. R. Luckhurst, Thin Solid Films 393, 40 (2001). 

12. R. Berardi, L. Muccioli, S. Orlandi, M. Ricci and C. Zan- 
noni, J. Phys.: Condens. Matter 20, 463101 (2008). 

13. L. A. Madsen, T. J. Dingemans, M. Nakata and E. T. 
Samulski, Phys. Rev. Lett. 92, 145505 (2004). 

14. B. R. Acharya, A. Primak and S. Kumar, Phys. Rev. Lett. 
92, 145506 (2004). 

15. A. G. Vanakaras, S. C. McGrother, G. Jackson and D. J. 
Photinos, Mol. Cryst. Liq. Cryst. 323, 199 (1998). 

16. R. van Roij, M. Dijkstra and R. Evans, Europhys. Lett. 
49, 350 (2000). 

17. R. van Roij, M. Dijkstra and R. Evans, J. Chem. Phys. 
113, 7689 (2000). 

18. M. Dijkstra, R. van Roij and R. Evans, Phys. Rev. E 63, 
051703 (2001). 



S. Varga et al.: Competition between capillarity, layering and biajciality in a confined liquid crystal 



13 



19. L. Harnau and S. Dietrich, Phys. Rev. E 66, 051702 (2002). 

20. T.J. Sluckin and A. Poniewierski, in Fluid Interfacial Phe- 
nomena, edited by C. A. Croxton (Wiley, New York, 1985). 

21. M. M. Tclo da Gama and P. Tarazona, Phys. Rev. A 41, 
11491152 (1989). 

22. T. J. Sluckin, A. Poniewierski, Mol. Cryst. Liq. Cryst. 179, 
349 (1990). 

23. Z. D. Zhang and W. J. Ye, Liq. Cryst. 36, 885 (2009). 

24. P. Shcng, Phys. Rev. Lett. 37, 1059 (1976). 

25. P. Shcng, Phys. Rev. A 26, 1610 (1982). 

26. T. J. Sluckin and A. Poniewierski, Phys. Rev. Lett. 55, 
2907 (1985). 

27. A. Poniewierski and T.J. Sluckin, Liq. Cryst. 2, 281 (1987). 

28. E. C. Gartland, P. Palffy-Muhoray and R. S. Varga, Mol. 
Cryst. Liq. Cryst. 199, 429 (1991). 

29. S. Kralj, S. Zumer, and D. W. AUender, Phys. Rev. A 43, 
2943 (1991). 

30. G. P. Crawford, D. W. AUender, and J. W. Doane, Phys. 
Rev. A 45, 8693 (1992). 

31. G. P. Crawford, R. Ondris- Crawford, S. Zumcr, and J. W. 
Doanc, Phys. Rev. Lett. 70, 1838 (1993). 

32. G. S. lannacchione, G. P. Crawford, S. Zumer, J. W. Doane 
and D. Finotello, Phys. Rev. Lett. 71, 2595 (1993). 

33. M. C. J. M. Visscnbcrg, S. StaUinga, and G. Vertogen 
Phys. Rev. E 55, 4367 (1997). 

34. P. Ziherl and S. Zumer, Phys. Rev. Lett. 78, 682 (1997). 

35. P. Ziherl, F. Karimi Pour Haddadan, R. Podgornik and S. 
Zumer, Phys. Rev. E 61, 5361 (2000). 

36. G. Panasyuk and D. W. AUender, J. Appl. Phys. 91, 9603 
(2002). 

37. P. Palffy-Muhoray, E. C. Gartland, and J. R. Kelly, Liq. 
Cryst. 16, 713 (1994). 

38. H. G. Galabova, N. Kothekar and D. W. AUender, Liq. 
Cryst. 23, 803 (1997). 

39. A. Sarlah and S. Zumer, Phys. Rev. E 60, 1821 (1999). 

40. C. Chiccoli, P. Pasini, A. Sarlah, C. Zannoni, and S. 
Zumer, Phys. Rev. E 67, 050703 (2003). 

41. G. Carbone, G. Lombardo, R. Barberi, I. Musevic and U. 
Tkalcc, Phys. Rev. Lett. 103, 167801 (2009). 

42. P. G. Dc Gcnncs, Langmuir 6, 1448 (1990). 

43. S. Kralj and T. J. Sluckin, Phys. Rev. E 50, 2940 (1994). 

44. M. Yamashita, J. Phys. Soc. Jpn. 72 1682 (2003). 

45. M. Yasen, M. Torikai and M. Yamashita, J. Phys. Soc. 
Jpn. 73 2453 (2004). 

46. M. Yasen, M. Torikai and M. Yamashita, Opto-Electron. 
Rev. 17, 89 (2009). 

47. M. Kio, M. Torikai and M. Yamashita, Opto-Electron. 
Rev. 17, 8 (2009). 

48. M. M. Telo da Gama, P. Tarazona, M. P. Allen and R. 
Evans, Mol. Phys. 71, 801-821 (1990). 

49. P. G. Ferreira and M. M. Telo da Gama, Physica A 179, 
179-198 (1991). 

50. D. Cleaver and M. P. Allen, Mol. Phys. 80, 253-276 (1993). 

51. C. Chiccoli, P. Pasini, F. Semeria, E. Bcrggren and C. 
Zannoni, Mol. Cryst. Liq. Cryst. 290, 237 (1996). 

52. C. Chiccoli, O. D. Lavrcntovich, P. Pasini and C. Zannoni, 
Phys. Rev. Lett. 79, 44014404 (1997). 

53. P. Pasini, C. Chiccoli and C. Zannoni, in Advances in the 
computer simulation of liquid crystals, P. Pasini and C. Zan- 
noni eds. (Kluwer, Dordrecht, 2000). 

54. G. D. WaU and D. J. Cleaver, Phys. Rev. E 56, 43064316 
(1997). 



55. Z. Bradac, S. Kralj and S. Zumer Phys. Rev. E 58, 
74477454 (1998). 

56. T. Gruhn, M. Schoen, J. Chem. Phys. 108, 9124 (1998). 

57. S. Varga and G. Jackson, Chem. Phys. Lett. 377, 6-12 

(2003) . 

58. H. Steuer, S. Hess and M. Schoen, Phys. Rev. E 69, 031708 
(2004) 

59. D. L. Cheung and F. Schmid, Chem. Phys. Lett. 418, 392- 
396 (2006). 

60. L. V. Mirantsev and E. G. Virga, Phys. Rev. E 76, 021703 
(2007). 

61. D. J. Cleaver and P. I. C. Teixeira, Chem. Phys. Lett.338, 
1 (2001). 

62. I. Rodn'guez-Ponce, J. M. Romero-Enrique and L. F. RuU 
Phys. Rev. E 64, 051704 (2001). 

63. I. Rodn'gucz-Poncc, J. M. Romero-Enrique and L. F. RuU, 
J. Chem. Phys. 122, 014903 (2005). 

64. D. de las Hcras, E. Velasco and L. Mederos, J. Chem. Phys. 
120, 4949 (2004). 

65. D. L. Cheung and F. Schmid, J. Chem. Phys. 120, 9185 

(2004) . 

66. H. Reich and M. Schmidt, J. Phys.: Condens. Matter 19, 
326103 (2007). 

67. P. I. C. Teixeira, F. Barmes, C. Anquetil-Deck and D. J. 
Cleaver, Phys. Rev. E 79, 011709 (2009). 

68. I. Rodn'guez-Ponce, J. M. Romero-Enrique, E. Velasco, L. 
Mederos and L. F. RuU, J. Phys.: Condens. Matter 12, A363 
(2000). 

69. A. Chrzanowska, P. L C. Teixeira, H. Ehrentraut and D. 
J. Cleaver, J. Phys.: Condens. Matter 13, 4715 (2001). 

70. P. I. C. Teixeira, A. Chrzanowska, G. D. Wall and D. J. 
Cleaver, Mol. Phys. 99, 889 (2001). 

71. D. dc las Hcras, E. Velasco and L. Mederos, Phys. Rev. E 
74, 011709 (2006). 

72. A. G. Vanakaras, M. A. Bates and D. J. Photinos, Phys. 
Chem. Chem. Phys. 5, 3700 (2003) 

73. Y. Martmez-Raton, S. Varga and E. Velasco, Phys. Rev. 
E 78, 031705 (2008). 

74. B. Mulder, Phys. Rev. A 35, 3095 (1987). 



